Bragg scattering of Cooper pairs in an ultra-cold Fermi gas 
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We present a theoretical treatment of Bragg scattering of a degenerate Fermi gas in the weakly 
interacting BCS regime. Our numerical calculations predict correlated scattering of Cooper pairs 
into a spherical shell in momentum space. The scattered shell of correlated atoms is centered at half 
the usual Bragg momentum transfer, and can be clearly distinguished from atoms scattered by the 
usual single-particle Bragg mechanism. We develop an analytic model that explains key features of 
the correlated-pair Bragg scattering, and determine the dependence of this scattering on the initial 
pair correlations in the gas. 
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Bragg scattering provides a high precision spectro- 
scopic technique that has been adapted from materials 
science to probe Bose-Einstein condensates [li yj. In 
condensate systems, signatures of soliton evolution [3|, 
phase fluctuations |j| , centre-of-mass motion Q , and vor- 
tex structure p, are accessible due to the velocity se- 
lectivity of Brag g spectroscopy. It has been proposed 
0, la, IE ll^ UM that Bragg spectroscopy of an ultra- 
cold Fermi gas can provide insight into the Cooper paired 
regime, and the transition through a Feshbach resonance 
to molecule formation. 

In this paper we develop Bragg spectroscopy as a probe 
of ultra-cold atoms by investigating Bragg scattering of 
a weakly attractive Cooper paired Fermi gas. Our cal- 
culations differ from existing theoretical treatments by: 
(i) providing explicit solutions for the time evolution of 
the matter field subjected to a moving optical grating, 
allowing direct observation of the dynamic response of 
the gas in momentum space, (ii) investigating in detail 
the large momentum transfer regime, where the atoms 
scatter well outside the Fermi sea, and, (iii) determin- 
ing the Bragg spectrum of the Fermi gas in an analo gou s 
way to the case for a Bose-Einstein condensate |3, ll^ , 
by calculating the momentum transfer per atom over a 
range of Bragg frequencies. The key result we report is 
Bragg scattering of correlated atom pairs via generation 
of a Bragg grating in the pair potential. 

Our theoretical treatment is based on a mean-field de- 
scription of a degenerate weakly attractive homogeneous 
Fermi gas. Two spin states are present in equal num- 
bers, with field operators '!/'-f(r, i) and ■!/'|(r,i), and the 
collisional interaction Hamiltonian is approximated by a 
number of single-particle terms. For convenience, the 
optical Bragg field is chosen so that it does not flip the 
particle spin. Implementing the Heisenberg equations of 
motion and the Bogoliubov transformation, 

k 

yields [l3, llJ, UM ^he time-dependent Bogoliubov de 



Gennes equations. 
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The Hartree potential is U{r,t) = V{^p}^{r,t)i(:a{'r,t)), 
the pair potential is A(r,t) = —V{il>-i{r,t)'ipi{r,t)), and 
Mk(r, t) and Vk{r, t) are the time-dependent quasi-particle 
amplitudes. We denote the Fermi energy by Ep — hujp — 
h^kp/2M, the atom mass by M, and the strength of the 
collisional interaction between fermions by V^ (V^ < 0). 
The Bragg field is T4pt(r, t) — Acos(q • r — ojt)/2, where 
the wave vector q is aligned with the x axis. 

The description of the collisional interaction in this 
problem requires a more realistic collisional potential 
than a contact potential, since we find that Bragg scat- 
tering is sensitive to the range of the potential. However, 
this sensitivity is weak for any realistic range, and we 
approximate the collisional potential using a contact po- 
tential with strength V and momentum space cutoff kc, 
where V and fcc are chosen to give the correct scattering 
length and correct momentum space range. While the 
static BCS properties are insensitive to kc |lflj, the total 
number of Bragg scattered atoms depends on the value 
of the cutoff fee, and in fact the number of pairs scattered 
vanishes on setting kc —> oo. Thus, kc must be chosen 
finite as above. Full details will be given elsewhere jlSi] . 

Details of our numerical method [l5j are beyond the 
scope of this paper. Briefly, the system ground state is 
found by solving iteratively for the self-consistent fields 
in the time-independent form of Eq. (|2Jl. The dynamic 
evolution of the gas is then determined by solving Eq. Q 
using a fourth order Runge-Kutta method [121 • We ex- 
plicitly evolve approximately 40, 000 k modes in order 
to describe the dynamics of a Bragg scattered three- 
dimensional zero temperature homogeneous Fermi gas. 
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FIG. 1: Column density, and pair correlation function, of a 
Bragg scattered three-dimensional homogeneous Fermi gas at 
zero temperature. The momentum distribution of the initial 
cloud is saturated on the chosen scale in order to observe the 
scattered atoms. Parameters are A = 1.80-Ef, q = 4.80fcF, 
t = 8.22/a;F, fee = SOfcp, fcpa = -0.427 [i.e., C/(0) = -0.256£;f 
and A(0) = 0M9Ef], and (a), (e) lu = 9.2a;F, (b), (f) lo = 
lO.Sojp, (c), (g) Lj = 14.0ajF, and (d), (h) lj = 28.4a;F. 



Figure n!^a)-(d) shows the monientuni space column 
density of a Bragg scattered Fermi gas for different Bragg 
frequencies to. The column density is J n{'k,t)dkz, with 
the number density at momentum Kk given by n(k, t) = 
C{(t>l{k,t)(j)a{k,t)), where 0a(k,i) = / ^^^(r, t) exp(-ik- 
r)(i'^r/(L)'^/^ is the field operator in momentum space, 
C is a constant chosen such that J n{'k,t)d^k = 1, and 
L^ is the computational volume. The two spin states 
are scattered identically because of our choice of a spin 
preserving Bragg field. 

The initial cloud, centered at k = 0, has an approxi- 
mate width given by the modified Fermi wave vector fcp, 
defined by h^k'^^ /2M = Ep - U{0) (e.g., [3). fn the 
presence of the Bragg field the atoms scatter by two dif- 
ferent mechanisms: (i) single-particle scattering, and (ii) 
correlated-pair scattering. 

In single-particle Bragg scattering an atom receives a 
momentum kick of nhq and energy nhiu, where n is an in- 



teger. A resonance condition selects primarily one value 
of n, and we consider u in the range of first order Bragg 
scattering (e.g., [l2|), where n = 1 is dominant. This re- 
sults in a cloud of atoms centered at k = q, as observed 
in Fig.nia)-(d). A faintly visible cloud at k = — q is also 
observed in Fig. ^a)-(d) due to off resonant scattering 
into the n = — 1 order. 

Figure ma)-(d) also shows scattering of atoms into a 
spherical shell centered at k = q/2. We refer to the 
mechanism responsible for the atom shell as correlated- 
pair Bragg scattering. Correlated-pair scattering has a 
frequency threshold denoted wthrcs, at which atoms are 
scattered to k = q/2 [see Fig. ^a)]. Above threshold 
[see Fig.^b)-(d)], the atoms are scattered into a spheri- 
cal shell and the shell radius increases with Bragg fre- 
quency. Atoms scattered into the shell come primar- 
ily from the Fermi surface, and are correlated about 
k = q/2, as demonstrated by the pair correlation func- 
tion CJq/2(k, i) = ((/)|(q/2-|-k,i)(/)|(ci/2 — k,t)), shown in 
Fig.^e)-(h) for the kz = plane [2l|. As well as the ring 
of correlation due to the scattered pairs, we also observe 
two rings of radius fcp centered at k = ±q/2, which rep- 
resent correlation between an atom on the Fermi surface 
of the initial cloud, and an atom in the first order Bragg 
scattered cloud. 

We determine the Bragg spectrum of the degenerate 
Fermi gas by calculating the momentum transfer per 
atom along the Bragg axis, V{t) = /[ftk-q]n(k, t)(Pk, for 
a range of Bragg frequencies. The Bragg spectrum at zero 
temperature is given in Fig. |2Ia), and is dominated by 
the broad single-particle resonance familiar from Bragg 
scattering of a Bose- Einstein condensate (e.g., |lS|). The 
single-particle Bragg resonance is due to two-photon scat- 
tering events that scatter atoms by the Bragg momentum 
hq. The resonance condition is well approximated using 
energy conservation arguments for non-interacting parti- 
cles, and by considering an atom scattered from momen- 
tum Kk^ to /i(kR -|- q) we obtain the resonance condition 



^ / 2 

2M ^^ 



2kR • q) 



(4) 



The single-particle resonance [see Fig. EJa)] is centered 
at w = hq'^/2M, and its width is Sui « 2hqk'p/M + A/h, 
where the first term accounts for the momentum width 
of the initial cloud, and the second term is due to power 
broadening (e.g., \\^). 

Correlated-pair Bragg scattering occurs on the red- 
detuned side of the single-particle resonance, leading to a 
shght asymmetry in the Bragg spectrum [see Fig. |2Ia)]. 
The correlated scattering, which gives rise to the dis- 
tinctive spherical shell of atoms in momentum space 
[see Fig. ^a)-(d)], depends critically on the presence 
of Cooper pairs, and can not be understood by the 
usual single-particle scattering mechanism. Correlated- 
pair scattering is associated with the formation of a mov- 
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tered pair is given by a; — Wthrcs = ^^rei/-^^' ^^ confirmed 
by our numerical calculations. 

We can understand some important features of Bragg 
scattering of correlated pairs with an analytic treatment. 
Due to the periodicity of the Bragg field, the solutions of 
Eq. ^ have the Bloch form, and can be expanded as 
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where n is an integer. The self-consistent potentials are 
periodic, with the translational symmetry of the Bragg 
field, i.e., 



C/(r,t) = E[/„(t)e™(i-'"-"*) 
A(r,t) = ^A„(t)e"('i-'"-"*). 



(7) 



FIG. 2: (a) Bragg spectrum of a three-dimensional homoge- 
neous Fermi gas at zero temperature. Parameters are A = 
1.80£f, q = 4.80fcF, t = 8.22/a;F, fcc = SOfcr, and (dashed) 
fcpa = and (solid) fcpa = —0.427. The non-interacting case 
is calculated with the initial momentum width of the gas mod- 
ified to agree with the interacting case, i.e., fcp = 1.12fcF. The 
vertical dotted line indicates lo = ojthros. (b) Relative change 
5n(i) = |A„(f)| — |A„(0)| of the pair potential Fourier com- 
ponents (+) n = \ and (o) n = [see Eq. Q] for the case 
fcpa = -0.427. 



ing grating in the pair potential, which is well approxi- 
mated by A(r,t) « Ao(i) -h Ai(i) exp[i(q ■ r - wi)]. The 
pair potential coefficients Ag and Ai, after the Bragg 
pulse, are shown in Fig. |2Ib), where we observe that the 
homogeneous term Aq is depleted in the region of the 
single-particle Bragg resonance, while the moving grat- 
ing amplitude Ai, is created over a slightly more ex- 
tended region. Cooper pairs (of zero centre-of-mass mo- 
mentum) scatter from the moving grating in the pair po- 
tential to centre-of-mass momentum ftq. At threshold, 
[see Fig. ^a)] the pair potential grating provides just 
sufficient energy to scatter each atom of a pair to a final 
momentum ?iq/2. Above threshold [see Fig.^b)-(d)], ex- 
cess energy provided by the pair potential grating is dis- 
tributed equally between the two atoms of a pair, and the 
individual momenta of the atoms become h{q/2 ± kj-ci). 
A resonance condition for correlated-pair Bragg scatter- 
ing can be obtained from energy conservation arguments 
to be 
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which for fcroi = gives the threshold frequency wtiucs — 
hq'^/4M - hk'p'^/M. Above threshold, the additional ki- 
netic energy h'^k'^^^/2M carried by each atom of a scat- 



Evolution equations for the coeffcients a^{t) and b^{t) 
can be derived from Eq. |2l to be 

zh^^ = tv^l{\^)al{t) + - {al^,{t) -I- ati(i)] (8) 



and 
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where hLu^-''i'k) = ± [/i^(k + nq)'^/2M - Ep] - nhuj. Ini- 
tially the only non-zero coefficients are Og (for |k| > k'p) 
and 6q (for |k| < k'p). The mean-field coefficients [/,„ 
and Am must be obtained self-consistently, in particular 



Amit) = -Vj2Y.''nitK-^it), 



(10) 



with V = r/(l — aT), where T is the low energy T-matrix 
and a = Mkj2TT^h^ [11 0. 

First order correlated-pair Bragg scattering is medi- 
ated by a moving grating in the pair potential, arising 
due to terms to = 1, n = 0, 1 in Eq. pU|l becoming non- 
zero. Those terms are seeded by single-particle Bragg 
scattering events in which an atom on the Fermi surface 
is scattered by momentum hq. It can be shown, from 
Eqs. © and ^, that single-particle transitions gener- 
ate coefficients a^ (for |k| > fcp) and b^i (for |k| < fcp). 
In the region of the Fermi surface, lk| « kp, Ai is gen- 
erated due to the formation of contributions a^^feg* ^^'^ 
Cq fo*!* . Physically this corresponds to seeding pair corre- 
lations about k = q/2. When an atom with momentum 



Skpi « SfcpicR, is scattered by a single-particle transition 
to fi(kR + q) , it remains correlated with its unscattered 
pair at momentum —Hk^. This reduces the system pair- 
ing about k = (Aq), and generates a pair correlation 
about k = q/2 (Ai) [seeFig.|2Ib)]. This correlation seed- 
ing can also be observed in 5q/2(k, f) [see Fig.nie)-(h)], 
where a point on the left most circle represents correla- 
tion between an atom on the Fermi surface of the initial 
cloud and its pair which has been scattered by q. 

Following the initiation of the pair potential grating, 
its subsequent development can be understood in terms 
of a truncated version of Eqs. (jSJ and Q, i.e., 
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which is appropriate for describing the scattered pairs. 
In_Eq. dnj, e;^/(k) = hwl'^'i'^) ± t/o, and Ai(t) = 
-V^^a^{t)b^^{t). The term a^h^^ only becomes sig- 
nificant if e = eg(k) — e!Li(k) « 0, i.e., when correlated 
scattering transitions conserve energy [uj > Wthrcs)- At 
threshold, the summands in Ai have a stationary phase, 
leading to enhancement of the grating amplitude Ai [see 
Fig-EIb)]. The thickness 5k of the spherical shell of scat- 
tered pairs can be estimated by assuming a frequency 
width r, determined by the Bragg pulse length (F w Tr/i), 
and setting 5e — hV ^ to find that 
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We have investigated the dependence of correlated- 
pair Bragg scattering on a range of system parameters. 
In Fig. n^b), ~ 0.2% of the atoms are scattered by 
correlated-pair Bragg scattering, and the number of scat- 
tered pairs grows linearly with the length of the Bragg 
pulse (until t ^ 70/wf)- Over the range —0.18 > k-pa > 
—0.69, the number of pairs scattered increases quadrat- 
ically with A(0) indicating the coherent nature of the 
scattering process. For kpa = —0.689 [all other param- 
eters as per Fig. ^b)] there are ~ 6% scattered pairs. 
However, we emphasize that for \k-pa\ ^ 1 the mean- 
field approach may not be quantitatively accurate (e.g., 
[llj,|20j). The number of correlated pairs scattered can be 
further increased by enhancing the single-particle scatter- 
ing processes that seed the pair potential grating, either 
by increasing the Bragg field strength A^ or by reduc- 
ing the Bragg wave vector q (to make the seeding more 
resonant). 

In conclusion, we have calculated solutions of the time- 
dependent Bogoliubov de Gennes equations for a zero 
temperature homogeneous three-dimensional Bragg scat- 
tered Fermi gas, in the regime where the momentum 
transfer is well outside the Fermi surface. We predict 
Bragg scattering of correlated atom pairs, which has 
a distinctive signature in momentum space, namely a 
spherical shell of atoms centered at half the usual Bragg 



momentum transfer. Correlated-pair Bragg scattering 
occurs via a Bragg grating formed in the pair poten- 
tial, and has a well defined frequency threshold on the 
red-detuned side of the familiar single-particle Bragg res- 
onance. We have developed an analytic model that ex- 
plains the mechanism by which the pair potential grating 
is generated, and observe that the number of scattered 
pairs is proportional to the square of the initial pairing 
field. 
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